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We study the NAy transition in the Dyson-Schwinger approach. The nucleon and A baryons are 
treated as quark-diquark bound states, where the ingredients of the electromagnetic transition cur- 
rent are computed self-consistently from the underlying dynamics in QCD. Although our approach 
does not include pion-cloud effects, we find that the electric and Coulomb quadrupole form-factor 
ratios Rem and Rsm show good agreement with experimental data. This implies that the deforma- 
tion from a spherical charge distribution inside both baryons can be traced back to the appearance 
of p waves in the nucleon and A bound-state amplitudes which are a consequence of Poincare co- 
variance. On the other hand, the dominant transition amplitude, i.e. the magnetic dipole transition 
form factor, underestimates the data by ~ 25% in the static limit whereas agreement is achieved at 
larger momentum transfer, which is consistent with missing pion-cloud contributions. We further- 
more find that the static properties of the form factors are not very sensitive to a variation of the 
current-quark mass. 
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I. INTRODUCTION 

Investigating the structure of the nucleon represents 
a challenging task in contemporary hadron physics from 
both experimental and theoretical perspectives. Unam- 
biguous testimony of the nucleon's complex non-pointlike 
structure is given by measurements of the electromag- 
netic transition into its lowest-lying resonance, the A 
baryon. The properties of the electromagnetic spin- 3 /2- 
spin-!/2 transition are experimentally accessible via pion 
photo- and electroproduction off nucleon targets. Preci- 
sion measurements of the N — > A7 transition over a wide 
range of Q 2 have become available only in the past decade 
at facilities such as JLab/CLAS, MAMI, and MIT-Bates. 
Experimental data exist now up to Q 2 ~ 8 GeV 2 ; for 
comprehensive reviews on the present experimental sta- 
tus we refer to Refs. [1, 2], 

While it contributes less than 1% to the total de- 
cay width of the Delta-isobar, the N — > A7 transition 
represents an indubitable test of the deformation from 
sphericity in both nucleon and A baryons. It is char- 
acterized by three transition form factors: the magnetic 
dipole form factor Ml, the electric quadrupole form fac- 
tor E2, and the Coulomb quadrupole form factor C2, 
which can be equivalently expressed in terms of helicity 
amplitudes. Although the magnetic dipole transition is 
dominant, an accurate extraction of the form factor ratios 
Rem = E2/M1 and R S m = C2/M1 reveals small and 
negative values. This constitutes a measure of the non- 
spherical distribution of the partons within the baryons 
involved in the reaction. 

The theoretical description of the A" — > A7 transition 
has been traditionally quite challenging. In a constituent - 
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quark model where three quarks move non-relativistically 
in an s wave, the transition magnetic moment is underes- 
timated by 25%, and the electromagnetic ratios are only 
non-zero if d waves are included [3, 4]. More sophis- 
ticated approaches have been developed over the years 
with varying degrees of improvement. The A~ — > A7 
transition has been studied in chiral effective field the- 
ory [5, 6], dynamical reaction models such as the Sato- 
Lee model [7, 8], light-cone QCD sum-rule analyses [9], 
large- Nc approaches [10, 11], the general parametriza- 
tion method [12], the cloudy bag model [13, 14], Skyrme 
models [15-17], vector-meson-dominance types of mod- 
els [18], relativistic quark models improved by chiral cor- 
rections [19], and other constituent quark models [20-23]. 
Systematic progress in recent years has also been made 
by calculations in lattice QCD [24, 25]. 

The nucleon and A deformation is tightly connected 
to the understanding of their internal structure. In 
Poincare-covariant approaches, the source of deformation 
from sphericity is attributed to the quark orbital angular 
momentum distribution among a hadron's constituents, 
predominantly arising from p waves in the hadron ampli- 
tudes, and to a much lesser extent from d waves [26-28] . 
In addition, the properties of the nucleon and its first res- 
onance are sensitive to chiral corrections and expected to 
be generated by an interplay between pure quark states 
and pionic clouds. 

In support of this, a QCD-motivated analysis of the 
quark-core contribution to the A~ — > A7 transition prop- 
erties is desirable. We employ the framework of Dyson- 
Schwinger equations (DSEs) and Bethe-Salpeter equa- 
tions (BSEs) which provides a well-established tool of 
investigation. In the quantum field theory of color- 
charged particles, dynamical chiral symmetry breaking, 
confinement and the formation of hadron bound states 
are phenomena which require a non-perturbative treat- 
ment. DSEs represent a fully self-consistent infinite set 



2 



of coupled integral equations for QCD's Green functions. 
The Dyson-Schwinger framework is particularly suitable 
to address these phenomena as it provides systematic ac- 
cess to both perturbative and non-perturbative regimes 
of QCD; see [29-31] for reviews. 

In the present approach, hadrons are studied via co- 
variant bound-state equations, see [32-34] and references 
therein. Meson properties emerge from solutions of the 
qq Bethe-Salpeter equation, while baryon properties are 
obtained from its three-body equivalent, the covariant 
Faddeev equation. Although of far bigger complexity, 
the Faddeev equation was recently solved for the nucleon 
and A— baryon masses [35-37] and utilized to calculate 
the nucleon's electromagnetic, axial and pseudoscalar 
form factors [28, 38]. By implementing a rainbow-ladder 
(RL) truncation, i.e. a dressed gluon-exchange kernel be- 
tween any two quarks, these investigations enable a direct 
comparison with corresponding meson studies. Meson- 
cloud effects in the chiral and low-momentum structure 
of hadrons are presently not accounted for, hence the 
framework aims at a description of the hadronic quark 
core. 

It is interesting that the nucleon and A masses and 
nucleon form factors obtained from three-body Faddeev 
calculations show little discrepancies to those obtained in 
a simplified version of the approach, the quark-diquark 
model. The latter represents an efficacious simplification 
of the Faddeev three-body problem to a Bethe-Salpeter 
two-body problem, where correlations in the qq scatter- 
ing matrix beyond the dominant scalar and axial-vector 
diquarks are neglected [26, 39]. The quark-diquark ap- 
proach is based on the observation that the attractive 
nature of quark-antiquark correlations in a color-singlct 
meson is also attractive for 3c quark-quark correlations 
within a color-singlet baryon. In this respect, we consider 
a first study of the N — > A7 transition properties within 
a consistent and well-established quark-diquark frame- 
work worthwhile and therefore adopt this simplification 
in the present paper. Moreover, this study completes 
previous investigations of quark-core contributions to the 
nucleon and A— baryon masses and form factors [40-44], 
where the backbone of the approach was entirely pro- 
vided by the quark-diquark model. 

The manuscript is organized as follows: in Section II 
we discuss the Poincare-covariant Faddeev approach to 
baryons and its quark-diquark simplified setup. In Sec- 
tion III we elaborate on the construction and properties 
of the electromagnetic transition current operator in the 
quark-diquark framework. In Section IV we present and 
comment on the results for the N — > A7 transition form 
factors and compare them to experimental data. Tech- 
nical details of our calculations are presented in Appen- 
dices A-D. We work in Euclidean momentum space and 
use the isospin-symmetric limit m u = ma- 



ll. QUARK-DIQUARK FRAMEWORK 

The description of the A^A7 transition properties in 
the Dyson-Schwinger approach requires knowledge of the 
nucleon and A bound-state amplitudes and their micro- 
scopic ingredients in terms of QCD's Green functions. 
The quantities which appear explicitly in the computa- 
tion of hadron wave functions and form factors are the 
dressed quark propagator and the irreducible qq, qq and 
qqq kernels. They encode the interactions at the quark- 
gluon level and, at least in principle, can be obtained 
through QCD's Dyson-Schwinger equations. 

In order to compute nucleon and A properties, these 
building blocks must be combined via covariant bound- 
state equations. In the present study we treat baryons as 
bound states of quarks and diquarks which amounts to 
a simplification of the three-quark problem. Irreducible 
three-quark interactions are neglected, and the qq inter- 
action is subsumed in effective diquark correlations. As a 
consequence, gluons appear only implicitly in quark and 
diquark propagators and quark-diquark vertex functions. 

In subsection II A we will briefly describe the basic in- 
put of the approach, namely, the rainbow-ladder ansatz 
for the qq kernel which amounts to an iterated gluon ex- 
change between two quarks. The subsequent calculation 
of diquark and baryon properties from that input is dis- 
cussed in II B and II D, respectively. 

A. Quarks 

The basic building blocks that enter the description 
of hadronic bound states are the dressed quark propaga- 
tor S(p) and the two-quark irreducible kernel JC. These 
are the ingredients that appear in a meson's Bethe- 
Salpeter equation and the covariant Faddeev equation 
for a baryon. The quark propagator satisfies a Dyson- 
Schwinger equation which can be solved numerically 
within a suitable truncation. On the other hand, the 
absence of a self-consistent solution for the quark four- 
point function necessitates an ansatz for the kernel K,. 

The construction of kernel ansatze is guided by vec- 
tor and axialvector Ward-Takahashi identities which en- 
sure electromagnetic current conservation as well as Gell- 
Mann-Oakes-Renner and Goldberger-Treiman relations 
at the hadron level [38, 45, 46]. These identities impose 
constraints on the structure of the kernel K, by relating it 
with the kernel of the quark DSE [47, 48]. The simplest 
ansatz to satisfy those constraints is the rainbow-ladder 
kernel: 

- 47ra(fc 2 ) 

AW/3/3' - ^2 p 1 k lota'IW: U) 

where T£ v = — k^k v /k 2 is a transverse projector with 
respect to the gluon momentum k, and Zi is the quark 
renormalization constant. Eq. (1) describes a dressed 
gluon exchange between quark and antiquark that re- 
tains only the vector part ~ 7^ of the quark-gluon vertex. 
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FIG. 1. (Color online) Quark DSE (3) in rainbow-ladder trun- 
cation. 



Its non-perturbative dressing, together with that of the 
gluon propagator, is absorbed by an effective coupling 
a(k 2 ) which is modeled. 

The second basic ingredient is the dressed quark prop- 
agator. It is expressed by two scalar functions, the quark 
wave-function renormalization 1/A(p 2 ) and the quark 
mass function M(p 2 ): 



S- 1 (p)^A(p 2 ) {ip + M(p 2 )) 



(2) 



The quark propagator satisfies the quark DSE which is 
illustrated in Fig. 1. Its interaction kernel includes the 
dressed gluon propagator as well as one bare and one 
dressed quark-gluon vertex. In rainbow-ladder trunca- 
tion that kernel becomes identical to Eq. (1) and the 
quark DSE reads: 

S aB (?) = Z 2 + m o) a B + / K-aa'P'f) S a 'p> (q) . (3) 



The bare current-quark mass mo appears here as the in- 
put of the equation and can be varied from the chiral 
limit up to the heavy-quark regime. 

The quark DSE exhibits dynamical chiral symmetry 
breaking if the kernel K. supplies sufficient interaction 
strength. Its consequence is a non-perturbative enhance- 
ment of the quark mass function M(p 2 ) in the low- 
momentum region that indicates the dynamical genera- 
tion of a constituent-quark mass scale. In principle, such 
strength would be generated through a self-consistent 
DSE solution for the gluon propagator and quark-gluon 
vertex that enter the quark DSE. In the rainbow-ladder 
truncation, that effect is provided by the effective cou- 
pling a(k 2 ) which appears in the kernel (1). The ansatz 
we choose is taken from Ref. [49] and plotted in Fig. 2: 



a(k 2 



TTT] 



A 2 



a vv (k 2 ) 



(4) 



The second term auv is only relevant at large gluon mo- 
menta where it dominates and is constrained by per- 
turbative QCD. The important part in view of hadron 
properties is the first term. It provides the necessary 
strength at small and intermediate momenta that trig- 
gers the transition from a current-quark to a dynami- 
cally generated constituent quark. It is characterized by 
two parameters: an infrared scale A that represents the 




k 2 [GeV 2 ] 

FIG. 2. (Color online) Effective coupling a(k 2 ) of Eq. (4), 
evaluated for A = 0.72 GeV and in the range rj £ [1.6, 2.0]. 



scale of dynamical chiral symmetry breaking, and a di- 
mensionless width parameter r\ that modifies the shape 
of the interaction in the infrared, cf. Fig. 2. 

In combination with the interaction of Eq. (4), the 
rainbow-ladder truncation has been quite successful in 
describing a variety of hadron properties. Upon set- 
ting the scale A via the experimental pion decay con- 
stant, its implementation in hadronic bound-state equa- 
tions describes pseudoscalar-meson, vector-meson, nu- 
cleon and A ground-state properties reasonably well, 
see [28, 37, 50-53] and references therein. In addition, 
their features show an overall insensitivity to a variation 
of the infrared properties of the coupling which is con- 
trolled by the parameter r\ [49, 52]. 

Progress in the light meson sector has also been 
made for axial- vector and pseudoscalar isosinglet mesons. 
Their properties are subject to substantial corrections be- 
yond rainbow-ladder which mainly come from the quark- 
gluon vertex [54-56]. Another important type of correc- 
tions are pion-cloud contributions which are important 
in the chiral and low-momentum structure of hadrons. 
They give rise to chiral singularities and non-analyticities 
associated with the opening of decay channels. Pion- 
cloud effects arc not implemented in the present approach 
which therefore aims at investigating the properties of the 
nucleon and A quark core. 

Eq. (4) remains the basic model input throughout this 
work. By setting the scale A and treating 77 as a parame- 
ter that reflects the model uncertainty, the properties of 
the effective interaction directly translate to the level of 
form factors. The quantities that appear in intermediate 
steps are computed self-consistently from this input so 
that no further model assumptions are necessary. 
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B. Diquarks 



The importance of diquark correlations in view of un- 
derstanding baryon structure has a longstanding his- 
tory [57, 58]. The strong attraction in the color anti- 
triplet diquark channel has motivated the interpreta- 
tion of baryons as bound states of quarks and scalar 
('good') diquarks. Such a picture reduces the rotational 
degrees of freedom in baryons and can explain the ab- 
sence of baryon resonances that are predicted by sym- 
metric quark models. Nevertheless, recent lattice results 
for excited baryon masses show more similarity with tra- 
ditional quark model predictions [59]. 

In that respect, one should be cautious in order not to 
oversimplify the problem. Diquark correlations certainly 
carry internal structure, and the inclusion of axialvector 
diquarks is crucial as well [26, 60]. In our setup, diquarks 
approximate the structure of the qq scattering matrix T 
that enters a baryon's Faddeev equation. It is expressed 
as the separable sum of scalar and axialvector-diquark 
contributions: 



T, 



aUfS'oc' ~ 2^ 
D 



(5) 



whose ingredients, in terms of diquark propagators D 11 " 
and quark-diquark vertices T^, are computed dynami- 
cally from QCD. Here, the label 'D' denotes the type of 
diquark. For scalar diquarks one has fi = v = and for 
axialvector diquarks fi, v — 1 ... 4. The total qq momen- 
tum is called P, and p, p' are the relative momenta be- 
tween the quarks. A bar on an amplitude denotes charge 
conjugation, cf. App. A. 

The implementation of the full T— matrix in a baryon's 
bound-state equation is equivalent to solving its covari- 
ant Faddeev equation without irreducible three-quark in- 
teractions, as implemented in Refs. [35, 37]. From that 
point of view, the quark-diquark model defined by Eq. (5) 
merely amounts to a truncation of the Faddeev equation. 
In fact, both approaches yield quite similar results for 
nucleon and A masses and nucleon electromagnetic form 
factors [28, 37, 53], which indicates that scalar and axi- 
alvector diquarks provide indeed the overwhelming con- 
tribution to the binding of these baryons. The question 
of baryon excitations, on the other hand, has not yet 
been addressed in full sophistication in either of these 
frameworks. A recent study of the Roper resonance in 
a simpler model suggests a quark-diquark structure as 
well [61]. 

While the dominance of qq correlations in the struc- 
ture of baryons may be a generic feature, the presence 
of explicit timelike diquark poles in the T-matrix is tied 
to the properties of the rainbow-ladder kernel (1). The 
rainbow-ladder truncation generates qq poles in the full 
T— matrix which are recovered by the ansatz (5) in terms 




FIG. 3. (Color online) The quark-diquark BSE, Eq. (9). 



of scalar and axialvector diquark propagators: 



D 00 (P 2 ) 



D^{P 2 ) 



1 



P 2 + ml 



1 P 



(6) 



p2 



where Tp V is again the transverse projector defined be- 
low Eq. (1). Colored diquarks disappear from the phys- 
ical spectrum beyond rainbow-ladder; hence, these poles 
are truncation artifacts [48]. Nevertheless, they provide 
access to a dynamical computation of onshell diquark 
amplitudes and their masses m sc , m av from the diquark 
BSE: 



= K. 



^olol' 'ft' 'ft 



[S(q + )T»( q ,P)S T (-q^] a ,p l , 



(7) 



with q± = q±P/2. Eq. (7) has a similar structure as the 
BSEs for pseudoscalar and vector mesons: upon tracing 
out the Dirac structure, they only differ by a color fac- 
tor. Details about the Lorentz-Dirac structure of the di- 
quark amplitudes and the solution method can be found 
in Refs. [34, 62]. 

Diquark correlations in a baryon are offshell. The off- 
shell structure of the qq scattering matrix in the sepa- 
rable approximation is primarily encoded in the diquark 
propagators. Inserting Eq. (5) in the scattering equation 
for the T— matrix yields an expression that allows for 
a consistent determination of the propagators. It reads 
schematically [34]: 



d-\p 2 ) = Tr / t d {g q -kt x )t 



D : 



(8) 



where Go is the disconnected product of two quark prop- 
agators and K, the rainbow-ladder kernel. It involves the 
onshell diquark amplitudes from Eq. (7) that must be 
equipped with appropriate offshell ansatze. The resulting 
diquark propagators reproduce Eq. (6) on the respective 
mass shells but deviate from free-particle propagators at 
offshell momenta. 
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C. Nucleon and A bound-state equations 

The simplification of the three-body problem to an ef- 
fective quark-diquark interpretation of baryons proceeds 
via omitting genuine three-quark interactions and em- 
ploying the separable pole expansion (5) in the quark- 
quark T-matrix. The resulting interaction between quark 
and diquark is an iterative quark exchange, where in ev- 
ery iteration step the spectator quark and one quark in- 
side the diquark exchange roles [26, 34, 39, 63]. 

The corresponding quark-diquark BSE is illustrated in 
Fig. 3 and reads 

rg(p,P) = J K^ DQ S(k q )D^'(k d )T^(k,P), (9) 

k 

where Tg are the quark-diquark amplitudes of the re- 
spective baryon, with B = N or A. The quark-diquark 
exchange kernel is given by 

K&q = r o( fc - k d ) S T (q) T«(p r , Pd ) . (10) 

The scalar and axialvector diquark amplitudes and prop- 
agators that appear in these expressions were explained 
in connection with Eq. (5). p q>d and k q>d are the external 
and internal quark and diquark momenta, p r and k r are 
the relative momenta that enter the diquark amplitudes, 
and P is here the total baryon momentum, cf. Fig. 3. 
The Dirac-Lorentz structure of the baryon amplitudes 
rg is analyzed in the following subsection. For details 
on the solution of the quark-diquark BSE we refer to 
Refs. [34, 64]. 

The treatment of light baryons such as the nucleon 
and A in the quark-diquark approach usually retains the 
lightest diquark degrees of freedom, i.e., scalar and axial- 
vector diquarks. The A baryon then involves only axial- 
vector diquark correlations whereas the nucleon contains 
both. As a consequence, the N A7 transition will not only 
involve axial-scalar diquark transitions but also axial- 
axial correlations. We will return to the implications of 
this feature in Section IV A. 



D. Nucleon and A amplitudes 

When expressed as two-body bound states of quarks 
and diquarks, the Dirac-Lorentz structure of the nucleon 
and A amplitudes Tg is considerably simpler compared 
to the full three-body approach of Refs. [35, 37]. Their 
derivation and partial-wave decomposition in the quark- 
diquark model has been outlined in Refs. [26, 65]. In 
view of exploring the impact of different quark orbital- 
angular momentum eigenstates on the properties of the 
NA'y transition, we briefly repeat the construction here 
for convenience. 

A two-body system is described by two independent 
momenta, the total momentum P and the relative mo- 
mentum p. Three Lorentz-invariant momentum variables 
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TABLE I. Orthonormal basis elements for the nucleon and A 
quark-diquark amplitudes that appear in Eqs. (11-12). 



can be constructed out of those: p 2 , z := p ■ P, and P 2 . 
The total momentum is onshell: P 2 = — M|, where the 
index B denotes the nucleon or A— baryon. The quark- 
diquark amplitudes can be expanded in a basis whose 
dimension is determined by Poincare covariance and par- 
ity invariance, together with the restrictions that follow if 
the baryon is taken on its mass shell. The respective de- 
composition for the scalar and axial-vector diquark parts 
of the nucleon amplitude reads: 

2 

T%(p 7 P) = £ /£V, z) r k (r, P) A + (P) , 

Y (11) 
r&(p, P) = ]T /£V, *) r£{r, P) 7 5 A + (P) , 

fc=3 

and for the A amplitude: 

k=l 

The Lorentz-invariant coefficients fj? depend on the two 
invariant variables p 2 and z and are solutions of the 
quark-diquark bound-state equations. We have factored 
out the nucleon and A spinors and work instead with the 
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TABLE II. (Color online) Rest-frame partial-wave decomposi- 
tion of nucleon and A— baryon in the quark-diquark approach. 
The basis elements are characterized by their scalar and ax- 
ialvector diquark content and their eigenvalues with respect 
to quark-diquark spin (s) and orbital angular momentum (I). 
The colored boxes highlight the dominant components; for 
example, the dressing functions associated with t£ and Tg 1 
are much smaller than the remaining p— wave contributions. 



positive-energy and Rarita-Schwinger projectors which 
satisfy 

A + (P)u(P,s)=u(P,s), 
V^(P)u u (P,s) = u^(P,s). [ ' 

They are given by 

A + = i(l + f), P^=A+(T^-i 7 ^7^), (14) 

where P = P/(iMb) is the normalized baryon momen- 
tum, Tp V = — p^pv i s a transverse projector with re- 
spect to P, and 7^ = T^^ v is the transverse 7— matrix. 
The projectors inherit the constraints from the spinors: 

fA + = A+, P^P^ = ^-p flv = 0. (15) 

Instead of p and P, the basis elements in Eqs. (11-12) 
can be equally well expressed through orthonormal mo- 
menta P M and r M := pr^, i.e., such that r 2 = P 2 = 1 
and r • P = 0. The dependence on the Lorentz invari- 
ants p 2 and z is then carried by the coefficients fj^ only. 
This simplifies the construction of an orthogonal basis 
and is also convenient for practical calculations, e.g. in 
the baryon's rest frame, where P and r are Euclidean 
unit vectors. 

The largest linearly independent set of basis elements 
for the bound-state amplitude T^, and is given in 
Eq. (Bl). On the baryon's mass shell, which is enforced 



by the properties (15) of the projectors, the following 
independent basis elements remain: 

rV {1,/}, 

T%: { 7 £, r», P"} x {1, f}, (16) 

: {6^, 7 £ r v , r"r u , P^r"} x {1, /}. 

These can be further orthonormalized and arranged ac- 
cording to their (quark-diquark) spin and orbital angular 
momentum content in the baryon's rest frame, cf. App. B. 
The resulting classification in s, p 7 d and / waves is illus- 
trated in Tables I and II. We emphasize that p-wave con- 
tributions to the bound-state amplitudes emerge quite 
naturally because of Poincare covariance. Those disap- 
pear in the non-relativistic limit [65] but have important 
consequences for the behavior of the form factors in Sec- 
tion IV. 

III. ELECTROMAGNETIC TRANSITION 

A. NA*y transition current 

We now turn to the general properties of the NAj 
transition current and its decomposition in terms of 
Lorentz-invariant form factors. The current can be gener- 
ically written as 

J^(P,Q) = P^(P J )i 75 r^(P,Q)A+(P i ), (17) 

where Pi and Pf are the incoming nucleon and outgoing 
A momenta, with P 2 = -M^ and Pj = -M\. They can 
be expressed through the photon momentum Q — Pf —Pi 
and the average momentum P = (Pi+Pf)/2. The onshell 
structure of the current is ensured by the projectors de- 
fined in Eq. (14), i.e., the positive-energy projector A + for 
the nucleon and the Rarita-Schwinger projector P pQ for 
the A-baryon. Eq. (17) is a matrix in spinor space; the 
usual current matrix element (Pf, Sf | J M | Pi, Si) is ob- 
tained upon contraction with the A and nucleon spinors 
from Eq. (13). The momentum dependence of the projec- 
tors implies that the 7— matrices contained in the Rarita- 
Schwinger projector P pa (P/) are now transverse with re- 
spect to Pf. We extracted an explicit factor 75 in Eq. (17) 
so that the remainder r QM , which will be specified below, 
has positive parity. 

Similarly to the nucleon and A bound-state ampli- 
tudes, the composition of the four-point function T atJ/ in 
Eq. (17) is determined by Poincare covariance. For its 
explicit construction it is again convenient to work with 
orthogonal momenta. This is not yet the case for P and 
Q because the non- vanishing N-A mass difference entails 
P • Q 7^ 0, cf. Eq. (C3). We take instead the component 
of P transverse to Q: 

P£ = T^P V = P" - (P • Q) , (18) 

and normalize it to unity: := Pt ■ Here, 

Tq 1 = 5^ - Q»Q V (19) 
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is the transverse projector with respect to Q. Together 
with the normalized photon momentum Q, the current 
is now characterized by two orthonormal four-momenta, 
K and Q (instead of P and Q, or Pj and Pf), which will 
simplify its structure considerably. 

Using this construction, the most general form of the 
vertex T al1 that is compatible with Poincare covariance, 
positive parity and current conservation can be written 
as (cf. App. C2): 



a//, 



(20) 



where 7^ is transverse to Q. It depends on three real and 
dimensionless form factors gi{Q 2 )- 

For comparison with experiment, it is more convenient 
to work with the Jones-Scadron form factors G* M {Q 2 ), 
G* E (Q 2 ) and G i ^{Q 2 ) which are related to the pion elec- 
troproduction multipole amplitudes at the A— resonance 
position and can be expressed in terms of helicity ampli- 
tudes [2, 66]. The respective decomposition of the vertex 



iuj 
2\Z 



(G* M -G* E ) l5 e a ^ s K^Q 



E 1 Q L K 



G* c Q a K^ 



(21) 



where we used the dimensionless variables 
Q 2 



2(M 2 A + M 2 N ) 



_ (M A ±M N ) 2 + Q 2 
A± - 2{Ml+M%) [22) 



as well as uj := ^/A+A_ and b := J § (1 + M^/Mn). 

We show in App. C2 that the vertices in (20) and (21) 
are equivalent when contracted with the projectors in the 
current matrix (17), and the relations between the gi and 
the Jones-Scadron form factors are stated in Eq. (C18). 

Eq. (21) is identical with the standard Jones-Scadron 
expression [2, 66] which is given in terms of the Lorentz 
structures 



£.afi~/8 pi pS 



M 2 A 



M 2 N 



iuj e 



K~<Q S 



^.aXyS pi pS gfiXpa pP pa 

{Mi + M 2 N y 

Q a (Q 2 PK-P-QQv) 



1 Q 1 K 



(23) 



2iujTQ a K f * 



These relations can be verified by expressing Pj and Pt 
through P and Q and subsequently in terms of the unit 
vectors K and Q via Eq. (C6). 



B. Electromagnetic current in the quark-diquark 
approach 

The computation of the A^A7 transition matrix of 
Eqs. (17) and (21) from its substructure in QCD re- 
quires a microscopic description of its ingredients. A 




FIG. 4. (Color online) General expression for the ./VA7 tran- 
sition current in the quark-diquark approach, see Eqs. (24) 
and App. D. 



systematic construction principle to derive the coupling 
of a hadron to an external current is the 'gauging of 
equations' method of Refs. [67-69]. The procedure was 
applied in [70] to derive the relevant diagrams in the 
quark-diquark system; recent discussions and applica- 
tions in the three-quark framework can be found in 
Refs. [28, 38, 71]. 

Applied to our case, the basic idea is that the AA7 
transition matrix element is the AA pole residue of the 
quark-diquark Green function that is struck by an exter- 
nal photon. If the current systematically couples to all 
internal constituents, which means that it has the formal 
properties of a derivative, electromagnetic current con- 
servation is automatically satisfied. The photon coupling 
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at the hadron level is thereby resolved to its interaction 
with quarks and diquarks. The resulting current has the 
generic form 



J M = r A G 



1 Q-DQ 



G r 



N ■ 



(24) 



which is detailed in App. D and illustrated in Fig. 4. 
and Ta denote the incoming and outgoing baryon ampli- 
tudes, Go is the disconnected product of a dressed quark 
and diquark propagator, and T M represents the photon's 
coupling to Go which yields the impulse-approximation 
diagrams (a-b) in Fig. 4. Those alone are however not 
sufficient: in order to comply with electromagnetic gauge 
invariance, one must take the coupling to the quark- 
diquark kernel into account as well. This leads to an 
exchange-quark term (c) and further diagrams that in- 
volve seagulls (d-e), i.e., the photon's coupling to the 
quark-diquark vertices. 

This decomposition has been frequently used in the 
computation of nucleon form factors by employing model 
ansatze for its ingredients [63, 72]. Building upon 
Refs. [39, 40], our goal in the present work is to compute 
these ingredients selfconsistently. We have outlined the 
calculation of quark and diquark propagators, diquark 
amplitudes, and nucleon and A bound-state amplitudes 
in the previous sections. In addition, we also compute 
the quark-photon and diquark-photon vertices that ap- 
pear in diagrams (a-c) , so that no further model input is 
needed for those objects. The only quantities that require 
a certain amount of modeling are the seagulls which we 
will discuss below. We note that the same method was 
recently applied for studying nucleon [73] and A form 
factors [27] in the quark-diquark approach and, upon ne- 
glecting the seagull terms, also the AiV7r transition form 
factor [44]. 

The central ingredient of Fig. 4 that describes the 
microscopic interaction with the photon is the dressed 
quark-photon vertex T^(k,Q). Ultimately, the diquark- 
photon and seagull terms would be resolved to the pho- 
ton's coupling with quarks as well. The properties of 
the quark-photon vertex will therefore be reflected in all 
form-factor contributions, either directly via diagrams 
(a) and (c) or implicitly in diagrams (b), (d) and (e). 
Here we compute the vertex self-consistently from its in- 
homogeneous Bethe-Salpeter equation [74]: 



a/3 



lC aa >w [S(k' + )T^k',Q)S(k'_)] a , , , 



(25) 



see App. D3 of Ref. [28] for details of the solution method. 
The equation features the same rainbow-ladder kernel of 
Eq. (1) that also appears in the quark DSE and diquark 
BSE. Electromagnetic gauge invariance is expressed by 
the vector Ward-Takahashi identity, 



which is satisfied for the resulting vertex. In combination 
with analyticity at Q 2 — > 0, it allows to write the vertex 
as the sum of a Ball-Chiu term [75], which is fixed by 
gauge invariance, and a further transverse contribution: 



L rp , 



(27) 



where A(k 2 ) and B(k 2 ) are the dressing functions of the 
inverse quark propagator S~ 1 (k) = il/t A(k 2 )+B(k 2 ), and 
we abbreviated 



F(k 2 )+F{k 2 _) 



A, 



F(k 2 + ) - F(k 2 _) 
k\ - k 2 _ • 



with k± = k ± Q /2. As a consequence of the term 
that appears in Eq. (25), the structure of the BSE implies 
the existence of a transverse vector-meson pole in the 
quark-photon vertex [74] , see [38] for a recent discussion. 
Indeed, the inhomogeneous BSE self-consistently gener- 
ates a p— meson pole in the vertex at Q 2 — —m 2 which 
dominates the timelike structure of the NA-f transition 
form factors. 

Having determined the quark-photon vertex, the 
diquark-photon vertices that appear in Fig. (4b) are com- 
puted selfconsistently as well. The diquark content of 
the nucleon and A amplitudes necessitates the inclusion 
of axial-axial as well as axial-scalar transition vertices. 
Applying the gauging procedure to the inverse diquark 
propagator of Eq. (8) resolves the diquark-photon ver- 
tices into their quark-gluon substructure, where the pho- 
ton couples once again to all constituents and thereby en- 
sures a Ward-Takahashi identity for the vertex [69, 76]. 
This leads to various quark-loop diagrams that involve 
the quark-photon vertex, seagull terms, and also the 
qq kernel. The construction is detailed in App. A. 7 of 
Ref. [34] and can be directly adopted here. The prop- 
erties of the quark-photon vertex and, in particular, the 
p— meson pole that appears in its transverse term, are 
thereby recovered in the diquark-photon vertices as well. 

Finally, current conservation in the quark-diquark ap- 
proach also requires the inclusion of seagull vertices M M,a 
in diagrams (4d)-(4e) . These terms describe the coupling 
of the photon to the offshell quark-diquark vertices. As 
they are another source of model uncertainty they war- 
rant further discussion. The seagull vertices M M ' Q satisfy 
a Ward-Takahashi identity similar to Eq. (26) which in- 
volves differences of scalar and axial-vector diquark am- 
plitudes [39, 70, 77]. In analogy to the Ball-Chiu con- 
struction for the quark-photon vertex, the seagull vertex 
can be written as the sum of a part that is fixed by the 
WTI and analyticity at Q 2 — > 0, augmented by a further 
unconstrained transverse piece: 



JW WTI 



(28) 



Q"im<?) = s-\k + ) - s-\k-) , 



(26) 



The first term Mfy^ is a rather lengthy expression, es- 
pecially once the full Dirac-Lorentz substructure of the 
diquarks is taken into account. Its detailed form is given 
in App. A. 8 of Ref. [34]; here we only repeat the generic 
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structure: 



M WTI = e_Mi + e+M 2 ~ e dq M 3 , 
M W ti = -e+Mi - e_M 2 + e dq M 3 . 



(29) 



The first and second lines correspond to the expressions 
that enter Figs. (4d) and (4e), respectively. The flavor- 
charge traces e± and e dq represent the photon coupling to 
the seagulls' quark and diquark legs and are determined 
in App. D 3. 

Since the properties of the quark-photon vertex in prin- 
ciple also dictate the structure of the seagulls, the trans- 
verse term in (28) must involve a p— meson pole as well. 
In absence of a selfconsistent solution for the seagulls, we 
use the following ansatz for the p— meson part: 



-g(x)T^M^ 



'WTI ' 



(30) 



where x = Q 2 /m 2 p . Purely longitudinal terms in (28) 
do not contribute to the ./VA7 current because of cur- 
rent conservation; hence, this simple form amounts to 
an overall factor (1 — g(x)) in the WTI-conserving part. 
The following phenomenological function was chosen in 
Ref. [40]: 



1 x 2 
9{x) = — 



-p 3 (l+x) 



g p 1 + x 



(31) 



where 



= \f2m p j f p is computed from the p— meson 
BSE. The parameter p 3 will be discussed below; it con- 
tributes to the overall model uncertainty that is, so far, 
induced by the width parameter r\ in the effective cou- 
pling (4) only. 
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FIG. 5. (Color online) Current-mass evolution of nucleon and 
A masses from their quark-diquark BSEs. The band corre- 
sponds to a variation of the width parameter r\ = 1.8 ± 0.2; 
the star denotes the experimental value. We compare to a 
selection of lattice data for nucleon [78-82] and A [78-80] . 



IV. RESULTS AND DISCUSSION 

We summarize the steps that we performed to obtain 
the results in this section. The basic model input is the 
effective quark-gluon interaction of Eq. (4) that enters 
the rainbow-ladder kernel of Eqs. (1). Its infrared part 
depends on the scale A = 0.72 GeV which we adjust to 
reproduce the experimental pion decay constant, whereas 
the width r\ — 1.8 ± 0.2 remains a parameter and mod- 
ifies the infrared shape of the coupling a(k 2 ). Using 
this input, we solve the Dyson-Schwinger equation (3) 
for the quark propagator, the diquark Bethe-Salpeter 
equation (7) for the scalar and axial vector diquark am- 
plitudes, and we compute the diquark propagators from 
Eq. (8). These quantities are subsequently implemented 
in the quark-diquark BSEs (9) from which we obtain the 
nucleon and A masses and bound-state amplitudes. 

The resulting masses for nucleon and A were obtained 
in Refs. [27, 40, 53] and are shown in Fig. 5 as a function 
of the squared pion mass. They are in reasonable agree- 
ment with experiment and lattice results. The results ob- 
tained in the quark-diquark model are in several respects 
similar to recent three-body calculations of nucleon and 
A properties [28, 35, 37]. The nucleon mass is practically 



identical in both approaches, whereas the A mass in the 
quark-diquark model is larger by < 5% and also shows a 
larger dependence on rj which is illustrated by the bands. 
Analogous observations also hold for the nucleon's elec- 
tromagnetic form factors which are quantitatively sim- 
ilar in both setups except that the quark-diquark ap- 
proach produces a larger model dependence. Note that 
the rainbow-ladder truncation does not dynamically gen- 
erate a A —J- Nir decay width, i.e., our result for the 
A— baryon describes a stable bound state. Associated 
non-analyticities which would appear for M/^—Mn > 
are therefore absent. 

Having established the ingredients of the electromag- 
netic transition current, we calculate its matrix elements 
from the diagrams in Fig. 4 which are discussed in detail 
in App. D. The current implements the selfconsistent so- 
lution for the quark-photon vertex (25) and the diquark- 
photon vertex, and for the seagull vertices we use the 
expression in Eq. (28). The seagulls, which are necessary 
for current conservation in the quark-diquark model, rep- 
resent another source of model uncertainty. Their trans- 
verse p— meson parts are modeled by Eqs. (30-31) and 
affect the larger-Q 2 behavior of the form factors. We take 
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FIG. 6. (Color online) Q 2 — evolution of the magnetic dipole 
form factor G* M in comparison with experimental data from 
Refs. [83-87]. The band denotes the model dependence as 
discussed in the text. 



a variation of the parameter p% € [0,0.15] into account; 
the central value of that interval was used in Ref. [40] 
to maximize agreement for nucleon electromagnetic form 
factors at larger Q 2 . The combined model dependence, 
stemming from the the seagull variation together with 
the rj dependence in the effective interaction, leads to 
the colored bands in Figs. (6-8). 

The Jones-Scadron form factors G* M (Q 2 ), G* E (Q 2 ) and 
Gq(Q 2 ) are finally extracted from the Dirac traces in 
Eq. (C19). Since the approach is Poincare-covariant, the 
results arc independent of the choice of reference frame. 
In order to avoid complex continuations for the radial 
momentum variables in the N and A bound-state am- 
plitudes, we work in the frame where the photon mo- 
mentum is purely real: Q = (0, 0, \Q\, 0) or, expressed in 
terms of the unit vectors defined in Section C, Q — e% and 
K — e^. The singularities in the quark and diquark prop- 
agators that enter the form factor integrals restrict the 
accessible domain of photon momenta to Q 2 < 2.5 GeV 2 , 
see App. D 2. This value is quite small and due to the 
quark-diquark description; a genuine three-body calcu- 
lation would be able to reach Q 2 values roughly twice 
as large. In addition, the kinematic dependence on the 
non-vanishing N-A mass difference also imposes a lower 
limit for Q 2 . In order to obtain results at Q 2 = 0, we ex- 
trapolate the form factor results at non-zero momentum 
transfer using Pade approximants. The extrapolation re- 
gions are indicated by the dashed margins in Figs. (6-7). 



A. Q dependence of the form factors 

The NAj transition current is determined by the 
three Jones-Scadron form factors G M (Q 2 ), G E (Q 2 ) and 





M N 


M A 


G* M (0) 


Rem{0) 


Rsm(0) 


Exp. 


0.94 


1.23 


3.02(3) 


-2.5(5) 




Calc. 


0.94(1) 


1.27(3) 


2.23(2) 


-2.3(3) 


-2.2(6) 



TABLE III. Results at the physical u/d mass compared to ex- 
periment. Nucleon and A masses are in units of GeV, G* M (G) 
is dimensionless, and the ratios Rem and Rsm are given in 
percent. The experimental values for G M (0) and Rem(0) are 
the PDG values [88]. The parentheses in our results indicate 
the combined model dependence as discussed in the text. 



Gq(Q 2 ) which are experimentally extracted from the 
multipole amplitudes in pion electroproduction [1, 2]. 
The process is dominated by a magnetic dipole transi- 
tion (Ml) which, in a quark-model picture, amounts to 
a spinflip of a quark and is encoded in the form fac- 
tor G* M {Q 2 ). Its static experimental value is G* M (Q) = 
3.02(3) [88]; experimental data exist in the range up to 
Q 2 ~ 8 GeV 2 . The remaining electric (E2) and Coulomb 
(C2) quadrupole contributions are much smaller and 
measure the deformation in the transition. They are ex- 
pressed by the form factors G* E (Q 2 ) and Gq(Q 2 ) which 
are usually related to the magnetic dipole form factor 
through the form factor ratios 



R 



EM 



Rsm 



\Q\ G* c 



2 -Ma G* 



(32) 



M 



where |Q| denotes the magnitude of the photon three- 
momentum in the A rest frame. It can be expressed in 
terms of Lorcntz-invariant variables via 



Igj = tJ 
2M A 1 + 25 



(33) 



where uj was defined below Eq. (22) and 5 is related to 
the N-A mass difference, cf. Eq. (C2): 



5 = 



Mj- 
2 (M 2 A 



Ml 



(34) 



Our result for the magnetic dipole form factor G* M (Q 2 ) 
is shown in Fig. 6. We find good agreement with exper- 
imental data above Q 2 ~ 1 GeV 2 , whereas the quark- 
diquark result underestimates these data by ~ 25% in 
the limit Q 2 = 0, cf. Table III. This is comparable 
to constituent-quark model predictions [19], where the 
long-standing discrepancy with the data has been at- 
tributed to missing meson-cloud contributions. Their 
impact has been studied with dynamical reaction mod- 
els [7, 92], where the 'bare' A resonance extracted from 
the Nj* — > Nir scattering amplitude corresponds to the 
quark-core contribution and meson-cloud effects are gen- 
erated via rescattering processes. In these analyses the 
pion cloud is sizeable and accounts for ~ 30% of G* M (0) . 
Similar conclusions have been found in the cloudy bag 
model [13, 14] or covariant chiral quark models [19]. 
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FIG. 7. (Color online) Form-factor ratios Rem(Q 2 ) and Rsm(Q 2 ) compared to experimental data from Refs. [85-87, 89-91]. 
The bands denote the model uncertainty as described in the text. 



The same interpretation applies to our framework as 
well. Pion-loop effects are not implemented in a rainbow- 
ladder truncation, and their inclusion would yield char- 
acteristic non-analytic structures in the chiral and low- 
momentum structure of form factors. Our result for 
G* M (Q 2 ) is consistent with various nucleon form-factor 
calculations in the Dyson-Schwinger/Faddeev approach, 
where similar discrepancies were interpreted as signals 
of missing pion-cloud effects [28, 73]. Typical exam- 
ples are the nucleon's charge radii which underestimate 
the experimental values but converge with lattice data 
at larger quark masses; or the low-Q 2 behavior of elec- 
tromagnetic form factors that shows missing structure, 
whereas one hnds reasonable agreement with experiment 
at larger Q 2 . A recent calculation of the nucleon's ax- 
ial charge finds such a discrepancy as well [38]. On the 
other hand, the nucleon's isoscalar anomalous magnetic 
moment n s , where leading-order chiral corrections can- 
cel, is accurately reproduced by the Faddeev calculation. 
These observations suggest to identify the rainbow-ladder 
truncated nucleon with the 'quark core' in chiral effective 
field theories. 

In fact, we observe analogous discrepancies in the low- 
Q 2 region for all three form factors G* M (Q 2 ), G* E (Q 2 ) 
and G* C {Q 2 ). Consequently, these deviations cancel in 
the ratios of Eq. (32), and Rem{Q 2 ) and Rsm{Q 2 ) show 
a good agreement with the experimental data which can 
be seen in Fig. 7. Considering the absence of pion-cloud 
effects in our description, this result is quite remarkable. 
Non-zero values for Rem(0) and Rsm(0) are usually at- 
tributed to the presence of quark orbital angular momen- 
tum in the NAj transition, either via ri-wave components 
in the nucleon and A wave functions, or caused by pion- 
cloud contributions. A large sensitivity to pionic effects 
has been found in coupled-channel analyses as well, where 
the bare ratios extracted from the electroproduction data 



are close to zero, so that Rem and Rsm would be almost 
entirely dominated by the meson cloud [7, 92]. 

A related question concerns the asymptotic behav- 
ior of the form factors Rem and Rsm and the scale 
where a perturbative-QCD description sets in. Dimen- 
sional counting rules and hadron helicity conservation 
predict [93] 

Rem 1 , Rsm const, for Q 2 —> 00 , (35) 

which is clearly not realized for the available experimen- 
tal data. Rem is negative and small and remains prac- 
tically constant in the entire Q 2 range whereas Rsm is 
negative but rises in magnitude at larger Q 2 . Taking into 
account the orbital motion of the partons in the nucleon 
and A wave functions leads to a double-logarithmic cor- 
rection in the ratio Rsm of Eq. (35) [94]. 

While the large-Q 2 region remains inaccessible to our 
analysis due to the kinematic singularity restrictions, 
we can investigate the impact of quark orbital angular- 
momentum correlations in the bound-state amplitudes. 
As a consequence of Poincare covariance, the nucleon 
and A wave functions do not only include s and d— wave 
components but also p waves and, in the case of the A, 
even / waves. The corresponding basis decomposition 
in the quark-diquark model, in terms of eigenvalues of 
quark-diquark spin and orbital angular momentum in 
the respective rest frames, is shown in Table II. While 
the interpretation of spin and orbital angular momen- 
tum can change in different frames, the decomposition 
itself is Lorentz-covariant, and basis elements which cor- 
respond to p waves in the rest frame carry one power of 
the relative quark-diquark momentum unit r M . Our fol- 
lowing use of the terminology 's waves' and 'p waves' will 
therefore refer to that basis decomposition. 

The analysis of the N and A dressing functions asso- 
ciated with the dimensionless basis elements in Table II 
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FIG. 8. (Color online) Decomposition of G* E (Q 2 ) and Rem(Q 2 ) according to the orbital angular-momentum content of nucleon 
and A amplitudes. The dashed bands are the results obtained with the s— wave basis elements in Table II and the dotted bands 
follow upon including the p— wave elements as well. The experimental data are the same as in Fig. 7. 



exhibits a clear hierarchy in their magnitude: s waves are 
dominant, p waves are suppressed, and d waves provide 
a tiny contribution. The same observation can be made 
in the three-body description of the nucleon amplitude, 
where s waves contribute roughly 2 /3 to the nucleon's nor- 
malization and p waves the remaining third, whereas the 
contribution from d waves is at the order of only one per- 
cent [28]. Thus, one would expect the dominant orbital 
angular-momentum effects in the form factors to come 
from p-wave contributions as well. 

In order to test the sensitivity of the NAj transition 
to different orbital angular-momentum correlations, we 
computed the Jones-Scadron form factors upon retaining 
s— wave elements only. This leads to a notably different 
behavior for the electric quadrupole transition G E (Q 2 ), 
shown in Fig. 8. Except for very small Q 2 , its result car- 
ries now a negative sign. The corresponding value for the 
ratio Rem starts off close to zero and rises almost linearly 
with increasing photon momentum. Although the upper 
Q 2 limit is certainly too small to allow any statements 
concerning the perturbative behavior, a trend towards 
the perturbative prediction Rem — >> 100% in Eq. (35) is 
visible. If p waves are included, Rem changes its behav- 
ior and becomes negative, thereby reinstating agreement 
with the experimental data. The inclusion of further d— 
and /—wave basis elements produces only a minor change 
in that result. This observation implies that the negative 
value for Rem is indeed an effect of quark orbital angular 
momentum; however, it predominantly owes to p— wave 
effects in the amplitudes which are already generated 
from the dynamics of the quark core as a consequence 
of Poincare covariance. This demonstrates the impor- 
tance of relativistic effects in the properties of the NAj 
transition. Such effects are missed in the non-relativistic 
quark model, or in a covariant description where p waves 



are not taken into account. 

On the other hand, we find no such behavior for the 
remaining form factors G* M (Q 2 ) and G^(Q 2 ). The result 
for G M (Q 2 ) obtained with s waves only is practically 
indistinguishable from Fig. 6 throughout the Q 2 range, 
which indicates that orbital effects do no play a role in 
the magnetic transition form factor. The same is true for 
the Coulomb ratio Rsm which exhibits a similar shape 
as the full result in Fig. 7, except for a broader model 
uncertainty. 

We finally comment on the importance of axialvector- 
diquark degrees of freedom in the nucleon amplitude. 
Due to its isospin- 3 /2 nature, the A-baryon consists exclu- 
sively of an isospin-1 axial vector diquark . If the nucleon 
were made of a scalar diquark only, the NAj transition 
would be a pure axial-scalar transition and the quark 
impulse-approximation in Fig. (4a) would not partici- 
pate. Table IV shows that there are sizeable contribu- 
tions coming from the interaction of the A with the ax- 
ialvector component in the nucleon. These correlations 
appear in all diagrams of Fig. 4 and contribute ~ 40% 
to G M (0). An even more pronounced effect is again vis- 
ible in G E (0) whose axial-scalar contribution is negative 
throughout the Q 2 range, whereas its axial-axial compo- 
nent is positive. The combination yields a positive value 
for G E (0) and thus a negative value for R E m(0)- 

B. Quark-mass dependence 

We now turn to the current-quark mass dependence of 
the NA<y transition form factors. It is accessible by mod- 
ifying the current-quark mass that enters the quark DSE, 
and by implementing the corresponding quark propaga- 
tor in all subsequent equations. The quark mass can be 
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FIG. 9. (Color online) Current-quark mass dependence of the three form factors G M , Rem and Rsm at vanishing momentum 
transfer. The dashed vertical lines indicate the position of the u/d and strange-quark mass, and the stars are the experimental 
values from the PDG [88]. As there is no ss pseudoscalar meson in nature, the value m sa = 0.69 GeV corresponds to a 
meson-BSE solution at a strange-quark mass m s = 150 MeV [46]. 





G M (0) 


GH0) 


Gc(0) 


A(ax) - AT(sc) 


1.27(4) 


-0.05(2) 


0.53(3) 


A (ax) - AT (ax) 


0.96(2) 


0.11(1) 


-0.11(3) 


Full 


2.23(2) 


0.06(1) 


0.42(6) 



TABLE IV. Extrapolated results for the Jones-Scadron form 
factors at vanishing momentum transfer, and their decom- 
position according to the diquark content of the nucleon's 
bound-state amplitude. The first line corresponds to the 
scalar-diquark components in the nucleon and the second line 
to its axial- vector contributions. 



results can provide insight in the properties of the quark- 
core contribution. Its insensitivity to a quark-mass vari- 
ation, together with the overall agreement of Rem and 
Rsm with experimental data at non-zero Q 2 , suggests 
that the bulk contribution to these ratios is generated 
by the dynamics of the quark core. Thus, the impact of 
genuine pion-cloud effects in the form-factor ratios might 
be smaller than expected and mainly appear close to the 
chiral limit. Naturally, this question demands further 
investigation. 



further mapped onto the pion mass by simultaneously 
solving the BSE for the pion. 

The results in Fig. 9 show the extrapolated values of 
G* M , Rem and Rsm at vanishing photon momentum. 
The error bands combine the accumulated uncertainties 
from the model parameter rj, the seagull variation, and 
the extrapolation error. We find that, at least within the 
model uncertainties, the form factors G* M (0) and Rsm{0) 
are mostly insensitive to the pion-mass evolution. The 
only quantity that shows a more pronounced quark-mass 
dependence is the ratio Rem (0) which slowly decreases in 
magnitude. This might be again due to p— wave effects 
in the nucleon and A amplitudes which diminish with 
increasing quark masses and mainly have an impact on 
the electric quadrupole form factor. 

Further information on the pion-mass dependent prop- 
erties of the form factors G* M , Rem and Rsm is available 
from chiral effective field theory [5, 6, 95], see [2] for a 
review. Depending on the treatment of the small scales 
772 and Ma — Mjv, those results show different charac- 
teristics and typically exhibit a sizeable dependence on 
the pion mass. A consistent observation is that Rsm 
logarithmically diverges in the chiral limit as a result of 
including pion loops. 

Since pion-cloud effects are missed in our approach, our 



V. CONCLUSIONS 

We investigated the N — > A7 transition properties in a 
Poincare-covariant quark-diquark approach. Since pion- 
cloud effects are not yet incorporated, our results describe 
the quark core contributions to the transition form fac- 
tors. They are obtained by employing a rainbow-ladder 
truncation at the level of the dressed quark propagator 
which amounts to a dressed-gluon exchange between the 
quarks inside the diquark. Thereafter, all quark and di- 
quark ingredients are self-consistently determined from 
the corresponding Dyson-Schwinger and Bethe-Salpeter 
equations and thus related to the fundamental quantities 
in QCD. 

We find several characteristic features with impact for 
future investigations of the N — > A7 transition. As a con- 
sequence of Poincare covariance, the electric quadrupole 
transition is dominated by p— wave admixtures to the 
nucleon and A bound-state amplitudes. Upon removing 
p waves, the E2/M1 ratio becomes positive and grows 
toward the perturbative prediction Rem — 1- The mag- 
netic dipole and Coulomb quadrupole transitions, on the 
other hand, are generated from s— wave contributions 
alone. The contribution from d waves is almost negli- 
gible in all cases. 

All form factors depend only weakly on the current- 
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quark mass. The dominant transition in the N — > A7 
reaction, the magnetic dipole Ml, could receive ~ 25% 
pion-cloud corrections in proximity to the chiral limit. 
On the other hand, the E2/M1 and C2/M1 ratios re- 
produce the experimental data quite well. These results 
disfavor arguments which promote the pion cloud and/or 
d— wave contributions as the prevailing missing effects for 
a proper description of the non-spherical character of the 
transition. Instead, this feature is already encoded in the 
nucleon and A quark core, partly via s— waves and par- 
tially via quark orbital angular momentum in terms of p 
waves. 

From a future perspective, in alignment with the lat- 
est state-of-the-art calculations for nucleon and A masses 
and form factors, our approach can be improved by 
abandoning the diquark ansatz in favor of a Poincare- 
covariant solution of the three-quark Faddeev equation. 
Apart from eliminating the model dependence on diquark 
properties, such an investigation would also be able to 
reach higher photon momenta than presently feasible. 
Furthermore, a consistently developed truncation beyond 
rainbow-ladder will methodically enable the inclusion of 
missing chiral corrections. This will elucidate the under- 
standing of nucleon resonances and the interplay between 
pion-cloud contributions and the core of quarks. 



where a are the three Pauli matrices. The charge con- 
jugation matrix is given by 



C 



4 2 
7 7 > 



C T = c f = c 



-c. 



(A3) 



and the charge conjugates for (pseudo-) scalar, (axial-) 
vector and tensor amplitudes are defined as 

T(p,P) :=CT(-p,~P) T C T , 
T a (p,P) := -Cr a (-p,-P) T C T , (A4) 
T 0a (p,P) :=CT al3 (-p,-P) T C T , 

where T denotes a Dirac transpose. Four-momenta are 
conveniently expressed through hyperspherical coordi- 
nates: 



( VT^^/i 
VT^Vi 



y 2 sine; 
y 2 cos< 



z 2 y 



\ 



(A5) 



and a four-momentum integration reads: 
00 1 



1 1 



(27T) 



42 I d P 2 P 2 



dz y/l 



1 2?r 

dy 
1 b 
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Appendix A: Euclidean conventions 

We work in Euclidean momentum space with the fol- 
lowing conventions: 

4 

p-q = ^2pkqk, p 2 =p-p, i> = p-~i. (Ai) 

A vector p is spacelike if p 2 > and timelike if p 2 < 0. 
The hermitian 7— matrices 7^ = (7^)^ satisfy the anti- 
commutation relations {7^,7^} = 2S IMl/ , and we define 

^ = --[7*, 7"], 7 5 = -7V7V ■ (A2) 
In the standard representation one has: 

7 * = ( — fe ) 7 4 = ( 1 ^ 7 5 = ( A 



Appendix B: Partial- wave decomposition of nucleon 
and A quark-diquark amplitudes 

This appendix serves as a supplement to Section II D. 
We will highlight the general structure of the basis ele- 
ments Tfc , t£ and t£ p that enter the quark-diquark ampli- 
tudes in Eqs. (11) and (12) as well as their partial-wave 
decomposition. 

To count the number of basis elements in the nucleon 
and A bound-state amplitudes, we observe that and 
are three-point functions with two spinor legs (for 
quark and nucleon) and a scalar or axial-vector diquark 
leg. is a four-point function with two spinor indices, 
an axialvector-diquark index fi and a further Rarita- 
Schwinger index v. The largest linearly independent set 
of basis elements for each case reads: 

T^: {\,f,f,ff}, 

Y%: K,r^}x{l,/,f,/|>}, (Bl) 

&>" x{t, f,r,fh 

K, r", P"} x {7?, r v , P v } x {1, f, f, /#>}. 

In case of the nucleon, the projector A+(P) in (11) ab- 
sorbs all occurrences of tf via Eq. (15) so that two in- 
dependent basis elements remain for the scalar quark- 
diquark amplitudes and six elements for the axial-vector 
contributions. The same applies to the A, where the 
Rarita-Schwinger projector in addition eliminates all in- 
stances of 7y and P" which leaves the eight basis ele- 
ments in Eq. (16). 
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To obtain the decomposition in Table I we need 
to construct quark-diquark spin and orbital angular- 
momentum operators. We will illustrate the procedure 
here for the case where only the dominant diquark am- 
plitudes 7 5 C and 7 A 'C are retained. While the total an- 
gular momentum j of the baryon, corresponding to the 
eigenvalues of the squared Pauli-Lubanski operator, is 
Poincare-invariant, the classification in terms of quark- 
diquark spin and orbital angular momentum is not. As 
exemplified in Ref. [28] in the three-body framework, one 
can retain a covariant notation by introducing formally 
covariant three-quark operators and L M which in the 
rest frame coincide with the spin and orbital angular- 
momentum operators. The square of then reads 

5 2 = |l«)l«)l + i {o T v ® o% ® 1 + perm.) . (B2) 

Applying (B2) to the basic scalar and axialvector di- 
quark structures 7 5 C and j^C yields 

SL W ' 7 y h 5 C) >y = [fl] aa , (l 5 C)^ , (B3) 

•S'aa' / 9/3'77'(7 Al C)/3'7' = 

3 (IT"! TM)! (YC)^. (B4) 

\ /J aa' 



and thereby defines effective quark-diquark spin opera- 
tors, namely the square brackets in the above equation, 
that act on the quark index a and the axialvector diquark 
index /i. They yield eigenvalues s(s + 1) for appropriate 
linear combinations of (16). In the scalar-diquark case 
(B3), only s = 1 /2 can appear. In the axial- vector di- 
quark case (B4), the spin operator is the sum of spin- 1 ^ 
and spin- 3 /2 projectors: 
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(B5) 



where A' 1 " = Tjf 



3 7t7t- 



Thus, for instance, basis 

elements that exclusively depend on the quantities P M or 
7y must carry s = 1 /2 as they are orthogonal to K^ v , cf. 
Eq. (15). 

The quark-diquark orbital angular momentum is en- 
coded in the operator 



involve up to three powers in the momenta r a . The cor- 
responding eigenstates of L 2 are given by 



1 = 1: r a 

1 = 2: 3r a rP-Tf 

1 = 3: 5rVV-S [a/37] (T£V) 



(B8) 



where <5[ q m denotes a symmetric permutation of the in- 
dices a, f3, 7. 

Reexpressing the relative-momentum dependence of 
the basis elements in terms of (B8) allows to read off their 
orbital angular momentum content. In the case of the 
nucleon these structures are already visible in Eq. (16): 

• 1, 7^ and P M carry no relative-momentum depen- 
dence and correspond to I = 0; 

• f, , 7^ f and P^f are proportional to r and there- 
fore p waves; 

• and r^f must be replaced by 3r M / — 7^ to obtain a 
d wave. 

The procedure is analogous for the A baryon, where 8^ v 
is an s wave, {7^ r v , P^r u , 8^ v f} are p waves and so on. 
Combined with the respective spin eigenstates, one ar- 
rives at the decomposition of Fig. II and Table I. The 
structure dictated by Eq. (B8) is not immediately ap- 
parent in the A covariants since we have exploited the 
properties (15) of the Rarita-Schwinger projector that 
appears in the full amplitudes. 



Appendix C: Electromagnetic transition current 

1. Kinematics 

The NAj transition matrix of Eq. (17) involves two 
independent momenta, the incoming nucleon momentum 
Pi and outgoing A momentum Pt. They can be ex- 
pressed by the photon momentum Q and the average 
momentum P: 



L 2 = 2 PT ■ d p 



( M- v 

[PtPt 



Pt T£ v ) dH dl 



(B6) 



7 



Pi, P 



P + P 



f 



(CI) 



with eigenvalues 1(1 + 1). In the baryon's rest frame, it 
reduces to the usual orbital angular-momentum operator 



L 2 = 2p ■ V p +p k {p ■ V p ) V£ - p 2 A p . 



(B7) 



Only the action of L 2 on the basis elements (r, P) is 
relevant since the dressing functions of Eqs. (11-12) are 
Lorentz- invariant, and L 2 applied to any of the Lorentz- 
invariants p 2 and z = p ■ P yields zero. Moreover, the 
eigenstates of L 2 can be determined independently of 
their Dirac structure since L 2 only acts upon the relative- 
momentum dependence. The basis elements in Eq. (16) 



with the inversion P l = P-Q/2 and P f = P+Q/2. Only 
Q 2 remains as an independent variable since the nucleon 
and A are onshell: P 2 = -M% and Pf = —M\. For the 
remaining Lorentz-invariant combinations one obtains: 



P z = - 



Ml + M 2 N Q 



PQ = - 



Ml -M 



N 



and Pi ■ Pf = P 2 — Q 2 /A. To simplify the notation, we 
abbreviate 



M 2 ;= Mj + Ml 



_ Ml-Ml 
' AM 2 ' 



AM 2 



(C2) 
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from which we obtain 

P 2 = -M 2 (1 + t) , 
PQ = -2M 2 S , (C3) 
Pi ■ P f = -M 2 (l + 2r). 

Expressed in terms of S and r, the quantities defined in 
Eq. (22) are given by 



A ± = 1±V \ — +t, w= V* 2 + i-(1 + t), (C4) 

and the normalized transverse average momentum X of 
Eq. (18) becomes 

^ = ^ = ^( P ' I+ 2^)' ™ 
Vice versa, P and Q depend on K and Q via 

P p = (icjif - , Q p = 2MVr . (C6) 

These expressions are ideally suited for evaluation in the 
the frame where Q — e% and -fT = are simply the 
Euclidean unit four- vectors. 



Structure of the current 



remove all occurrences of 7" or P a in (C7) as well. One 
remains with the following four basis elements: 



r 



Q a P p , a "*, Q a Q p 



(Gil) 



each of which carries a Lorentz- invariant form factor that 
depends on Q 2 . Current conservation Q^J P ' P = leads 
to a linear relation between them and implies that J p - p 
can only involve basis elements that are transverse to Q 
in the index fi. The same effect can be achieved by ap- 
plying a transverse projector with respect to the photon 
momentum to (Cll). With an additional normalization 
of the momenta, this results in the expression given in 
Eq. (20): 



Tiau 
■ Q 



(C12) 



Next, we want to verify the equivalence of the two 
forms for the transition current obtained from Eq. (20) 
with Eq. (21). To this end, we must show that the ten- 
sor structures that appear in (21) can be reexpressed by 
those in (20) when contracted with the Rarita-Schwinger 
and positive-energy projectors. For notational conve- 
nience we abbreviate F pa = P pQ (P/) and A; = A + (P l ) 
in the following. 

The structure Q a _ftT M attached to requires no fur- 
ther examination as it appears in both equations. The 
relation 



The generic structure of the matrix-valued N A7 tran- 
sition current is given in Eq. (17). It involves T p ' a (P, Q) 
which is a four-point function of positive parity, with 
two spinor and two vector indices (jj, is the photon and a 
the Rarita-Schwinger index). Its most general Poincare- 
covariant tensor structure involves 40 elements: 



ryOl p/j pa pp. 

^ a Q p P a Q^ 

flap pay* Q*pv. 

Q a j p Q a Q p 



x{l,f, % [?,$]}■ (C7) 



Only a few of those will survive on the N and A mass 
shell which is enforced through the projectors in Eq. (17). 
Using the properties 



A + (P f )f f = iM A A + (P f ), 
f i K + {P i ) = iM N K + {P i ) 



(C8) 



of Eq. (15) allows to systematically eliminate the struc- 
tures f, $ and [f, so that in the second bracket of 
Eq. (C7) only the unit matrix remains as an indepen- 
dent element. Moreover, the properties of the Rarita- 
Schwinger projector, 



V pa (P f ) 7 Q = F pa (P f ) P? = 0, 



and thus 



V pa (P f )P a = -V pa (P f ) 



Q a 



(C9) 



(CIO) 



F pa K a = --{8-t) F pa Q a 



(C13) 



that follows from (C5) and (CIO) allows to relate the 

(C14) 



tensorial structure of G%, 



rpCtJ rpffl _ rpCtfJ, j^a Tf fl 



to Tq^ and Q a K p as well. The remaining term including 
the Levi-Civita symbol can be written as 



+ K p 1 ^-K a 1 p $ 



(C15) 



cf. Eq. (A14) of Ref. [28]. The terms with 7^ amount to 

Vflsl% = -V p f a lsQ a $, (C16) 

via Eq. (C9), and those featuring K a can be again re- 
duced using (C13). By expressing K and Q through Pj 
and Pf and using the relations (15) one further derives 



A/75 It 



A; 



(C17) 



A/75 



A, : . 
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Putting pieces together, one finally obtains: 
F pa e a ^ 5 K~<Q s A, = 



f 75 



LO 



T" p + Q a K^ 



28 



A, 



This expression features the same tensor structures as 
Eq. (20), and therefore the current constructed from 
Eq. (21) is indeed a reparametrization of Eq. (20). 

Taking into account the full set of terms in Eq. (21) 
finally allows to read off the relation between the §i{Q 2 ) 
from Eq. (20) and the Jones-Scadron form factors G* M , 



G* E and G* c : 



9i 



2Ai 



(G\ 



M 



92 = b 



53 



2A_ 



(G*m ~ G E ) — 



G*e) > 



c 



{5-t)G% 



G*e 



(C18) 



where b 



(1 + M A /M N ). 



3. Extraction of the form factors 

The final issue to address is how the form factors are 
extracted from the A^ A7 transition matrix (17) once J p ' p 
has been computed from its underlying dynamics. The 
simplest Lorentz scalars that involve J p ' p are obtained 
from the following Dirac traces and momentum contrac- 
tions: 

Tr[ 75 J^- 



7 



Tr [75 J p ' p ] T pp 



(C19) 



si = 

S2 = 

s 3 = Tr [75 J p ' p ] K p K p , 
and their relation with the Jones-Scadron form factors is 

3V1 - 4<5 2 \\+ V^VT+~2£ 



G\ 



M 



G% 



VI - 4£ 2 
Abuj 



S2 



S2 



5- 



si 



tVI + 25 



5- 



■ si 



(C20) 



3V1 - 45 2 A+ (1 + 25) 



S3 ■ 



Aibuj 2 5 — t 
Eqs. (C19) are most conveniently evaluated in the 



frame where 



e% and 



Here the [i = 3 com- 



ponent of J p,p is longitudinal to the photon momentum 
and vanishes because of current conservation (as long as 
the current is conserved microscopically). Therefore, the 
traces become 

51 = Tr[ 75 >V], 

52 = Tr [ 75 J 1 ' 1 ] , (C21) 
.S3 = Tr[ 7 5</ 4 ' 4 ], 

where i is summed over i = 1,2. 



Appendix D: Transition current in the 
quark-diquark framework 

1. Current diagrams 

In the following we collect the ingredients of the A^A7 
transition current matrix in the quark-diquark model 
which are depicted in Fig. 4. The explicit form of the 
current is given by a sum of impulse-approximation di- 
agrams (upper two panels in Fig. 4) and two-loop con- 
tributions which represent the photon's coupling to the 
quark-diquark kernel (lower three panels): 

j p - p = J f p «(p f ,p f )(x q + x dq r^re N (p l ,p l ) + 
f^.P/jx&^rftfo.p,)- (di) 



The quark-diquark amplitudes and are the so- 
lutions of the nucleon and A bound-state equations (9). 
Pi and Pf are incoming and outgoing on-shell momenta 
and Q = Pf — Pi is the photon momentum, cf. Eq. (CI). 
The relative momenta pi and pf are independent loop 
momenta in the two-loop diagrams; in the one-loop dia- 
grams they are related to each other: pf —pi = (1 — £) Q 
for the quark diagram and pf —pi — — £ Q for the diquark 
diagram, where £ € [0, 1] is an arbitrary momentum- 
partitioning parameter which must be specified prior to 
solving the quark-diquark BSEs. a, j3 — 1 ... 4 are the 
Lorentz indices of the axialvector diquark in the A and 
nucleon amplitudes, and (3 = corresponds to the addi- 
tional scalar-diquark component in the nucleon. 
The ingredients of Eq. (Dl) are given by 

= S( P+ ) T p (p, Q) S(p-) D^(k) , (D2) 
= S(p) D aa ' (k + ) Ttg? (k, Q) , 

x gct0 = D aa' \k+)S{p+)K fl ' a ' ' ' S(p-)DW{k-) 

and depend on the quark-photon vertex r p and the 

diquark-photon vertex r p 3 '" /3 , where the latter includes 
an axial-axial contribution and an axial-scalar transition 
component for j3 = 0. The quark and diquark momenta 
are: 



p- =Pi + £Pi, 
p+ =Pf + £ p f 



k- = -Pi + (1 - Pi , 

k+ = - Pf + (i-OPf , 



and we denote the average quark and diquark momenta 
that appear in the vertices by p = (p + + p_ ) /2 and k = 
{k + + k_)/2. 

The 'gauged' kernel AT p > a ^ contains the exchange- 
quark diagram and the seagull vertices M^' a : 



K p > af) = (K EX + K SG + K ; 



(D3) 
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with 



complex plane, 



K 



At, a/3 
EX ■ 



r^(r + ,fc_) S(q + )mq,Q)S(q-) T a (r_,k + ) 



K^=M^(r' + , fc_,Q) S(g+) T F°(r_, fc+) , 

^ = r ^ (r+ j fc _ } s{q _ )T m ,a {r ,_ > fc+ 5 Q) ^ (D4) 



and momenta: 



q± = k±- p T 



r± = 



P± - Qt 



p± - q± 



The color-flavor traces of the various contributions in 
Eqs. (D2-D4) are detailed in App. D3. 



2. Kinematic restrictions 

The integrand of the iVA7 transition matrix (Dl) 
contains dressed quark and diquark propagators which 
are computed selfconsistently from Eqs. (3) and (8). 
In rainbow-ladder truncation, the quark propagator ex- 
hibits singularities in the timelike complex plane and 
the diquark propagators have poles on the timclikc real 
axis. Without taking residues into account explicitly, 
that singularity structure restricts the accessible momen- 
tum phase space in the integrals (Dl) and leads to kine- 
matic limitations in the Q 2 range of the resulting form 
factors. 

The problem can be accessed by separating the real 
and imaginary parts of the quark and diquark momenta 
that enter the diagrams. As we have noted in the discus- 
sion of Eqs. (C6) and (C21), it is advantageous to work 
in the frame where the normalized momenta Q and K 
that describe the current are the Euclidean unit vectors 
Q = ez and K — e±. For spacelike values of Q 2 , the 
photon momentum is then purely real and the trans- 
verse average momentum P? is imaginary. This frame 
also suggests itself via the impulse-approximation rela- 
tions Pf — pi ~ Q from the previous subsection: if the 
photon momentum is real, both relative momenta pi and 
Pf can be chosen real as well (the loop momentum is al- 
ways real) and it is sufficient to determine the N and A 
amplitudes for pf, p 2 e R + . 

The quark and diquark momenta depend on the mo- 
mentum partitioning parameter £ £ [0, 1] which can 
be used for maximizing the available momentum phase 
space. We exemplify the procedure for the quark mo- 
mentum jj_ = pi + £ Pi. In the frame where the photon 
momentum is real: pi G K 4 , and the only imaginary con- 
tribution in Pi = P — Q/2 comes from the transverse 
component Pt, cf. Eq. (C6). Thus, p- can be written as 



P- 



(D5) 



where i? M is purely real and irrelevant for the further 
discussion. p_ describes the interior of a parabola in the 



t±i 



t e E+ 



(D6) 



that intersects the real axis at (p 2 _) max = — (M £ui) 2 / t . 
M is here the average N-A mass defined in Eq. (C2). 
On the other hand, the nearest singularities in the quark 
propagator (whether timelike-real or in the timelike com- 
plex plane) can be used to construct another parabola 



(t ± im q ) , t £ R+ 



(D7) 



which defines the quark 'pole mass' m q . Hence, in order 
to sample the quark propagator only in the kinematically 
safe region, one has to obey the condition 



> 



( mq V 



(D8) 



where we have used Eq. (C4) in the second step. Because 
of the non- vanishing N-A mass difference, expressed by 
5 0, this puts constraints on r = Q 2 /(AM 2 ) from both 
above and below. 

The analysis can be repeated for any quark or diquark 
momentum that enters the integrands (D2) and yields 
the combined relation 



+ T <h{0 2 -l, 



where 

h(^) = min 



m d 



CM' \1-2£\M' (l-OM 



(D9) 



(D10) 



with = min(77i, sc , m ax , 2m q ). In view of maximizing 
the kinematical region in Eq. (D9), we are interested in 
the value £o that maximizes h(£). For m q < < 2m q , 
that value is given by 



£ = (l + m d /m 9 )- 1 e[l/3...1/2], 



(Dll) 



which leads to h — h(£ ) — (ni q +m d )/M. ViaEq. (D9), 
the resulting lower and upper limits for Q 2 become 



r± 



hl-l 



1±W1- 
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i 



(D12) 



Inserting the numerical values for the masses m q , m sc , 
m ax , Mjv and yields a kinematically safe region that 
is centered around Q 2 ~ 1...1.3 GeV 2 , depending on 
the model parameter 77 in the effective coupling. 

We conclude with two remarks. First, while it would be 
possible to choose two different momentum partitioning 
parameters £at and £a in the calculation of the nucleon 
and A amplitudes and form factors, doing so would not 
relax the constraint in Eq. (D10): the optimal choice cor- 
responds to £/v = £a- The second remark concerns the 
choice of reference frame. The current diagrams could be 
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equally evaluated in the A rest frame where P/ = e± and 
where p 2 can be chosen as the real loop momentum. In 
that case, one has access to the limit Q 2 — and also to 
the timclikc region whereas the largest accessible space- 
like value is smaller, Q 2 ~ 1 GeV 2 . On the other hand, 
Q and Pi will become complex which necessitates a com- 
plex continuation of the nucleon bound-state amplitude. 
Such a procedure has been applied to compute the NAtt 
transition in the current framework [44]. 



3. Color and flavor traces 

In this appendix we collect the color and flavor traces 
that appear in the quark DSE, the various bound-state 
equations and the NAj transition matrix element. 

Each quark-gluon vertex in the kernel (1) is equipped 
with a color factor Ai/2, where the Aj are the eight Gell- 
Mann matrices. The color factors for the diquark am- 
plitudes are given by Sabc/VG and those for the N and 
A quark-diquark amplitudes by Sac I v3, where A, B are 
quark indices and C is the diquark index. Tracing the 
color structure leads to a prefactor 4 /3 in front of the inte- 
gral that appears in the quark DSE (3), a factor — 4 /3 for 
those in the meson BSE and the inhomogeneous BSE (25) 
for the quark-photon vertex, and a prefactor 2 /z in the 
diquark BSE (7). The color traces in the form-factor di- 
agrams of Fig. 4 yield +1 for the impulse-approximation 
diagrams and —1 for the exchange and seagull diagrams. 

Considering the flavor case, we work in the SU(2) 
isospin-symmetric limit with two degenerate quark fla- 
vors u = ( J ) and d = ( 5 ) . Out of those one can con- 
struct the following diquark flavor matrices: 



( i&2 



2 ' V2' 2 



(D13) 



where So corresponds to the isoscalar (scalar) diquark and 
the three Sj, with i = 1, 2,3, to the isotriplet (axialvector) 
diquark. The Oi are the Pauli matrices, and the flavor 
matrices are normalized to Tr{sJ Sj} = Sij. 

The flavor factors for baryons in the quark-diquark pic- 
ture are given by the Clebsch-Gordan coefficients accord- 
ing to the respective diquark content of the baryon. For 
proton and neutron one obtains 



P = 
n = I d 



|d,-,/±u,0 



0, Jfd, - A /|u 



(D14) 



where the first entry represents the isoscalar diquark and 
the remaining three the contributions from the isovec- 
tor channel, in the same order as Eq. (D13). For exam- 
ple, the two mixed-antisymmetric and mixed-symmetric 
flavor tensors of the proton's three-quark amplitude are 
s <8> Po and J2i s » ® Pi: respectively. On the other hand, 
the A-baryons do not contain any contribution from the 



scalar diquark, and the corresponding Clebsch-Gordan 
construction for the A + and the A (only these two are 
relevant for the NAj transition) yields 



A+ = 
A°=( i). 



§d. 



u, 



Id: 



(D15) 



The flavor traces in the nucleon and A quark-diquark 
BSEs (9) depend on the type of involved diquarks. This 
leads to combined color-flavor factors c' 00 * 1 = — c^ a ^ — 
— 1/2 and c^ ") = = v^/2 for the nucleon, whereas 

in the case of the A the color-flavor trace is given by 

c (a« = _l. 

The flavor factors for the form- factor diagrams of Fig. 4 
are obtained by combining Eqs. (D13-D15) with the 
quark charge matrix Q = diag(q Ul qd) which is attached 
to the quark-photon vertex. In the case of the pj —> A + 
transition, the three contributions in Fig. (4a-c) yield the 
following traces: 



]T(A+)t Q Pi , J2 L>+) f Pi] 2 Tr i s I s o Q f }, 

i ij 

^(A+)t Sj Qtslp,. 



(D16) 



Only the axial- vector diquark contributes to the A (hence 
i = 1,2,3) whereas the nucleon has both scalar (j = 0) 
and axial- vector diquark components (j = 1,2,3). 

The charge-flavor structure of the seagull diagrams in 
Fig. (4d-e) appears in combination with their Dirac com- 
ponents via Eq. (29). The flavor factors read: 



e_=5>+)t s ,Qtslp,, 

ij 

£(At)tQs iS | Pi> 



&dq = 



(D17) 



£[(A+)t Sj si Pj ]2Tr{s]s,Qt} 



for the seagull vertices in Fig. (4d), and 
e + = ^(A+)ts 3 s!Q Pj , 

ij 

e_=^(A+)ts,Qtst Pj , 



(D18) 



Pj 



2Tr{ S T s, Qt} 



for the conjugated seagulls that appear in Fig. (4e). 

In combination with the color traces +1 for the 
impulse- approximation diagrams and —1 for the ex- 
change diagrams, the final result for the p"f — > A + tran- 
sition matrix element becomes 



A 



V2 
3 



[Jq + JdQ + ^EX + ^SG + <%g] : 
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where the individual contributions are given by In this notation, the A, denote the Dirac parts obtained 

from Eqs. (D2) and (D4), and the superscripts S and A 
= > refer to the scalar or axial-vector diquark content in the 

_ _^aa _|_ y/^y^AS (D19) outgoing A (left) and incoming nucleon (right) ampli- 

i • ; dq dq tudes. The labels {1,2,3} for the seagull contributions 

>^ex = | Agx — 2 ^ex > correspond to those in Eq. (29). Identical traces follow 



2 EX 2 

±(A AA + A^ - 



j SG = i (A aa + a aa } _ a aa _ ^3 (a as _ A A S) for the process n 7 -> A . 
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